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Abstract — This paper deals with physical layer network-coding 
for the four-way wireless relaying scenario where four nodes A, 
B, C and D wish to communicate their messages to all the other 
nodes with the help of the relay node R. The scheme given in the 
paper is based on the denoise-and-forward scheme proposed first 
by Popovski et al. in |1|. Intending to minimize the number of 
channel uses, the protocol employs two phases: Multiple Access 
(MA) phase and Broadcast (BC) phase with each phase utilizing 
one channel use. This paper does the equivalent for the four-way 
relaying scenario as was done for the two-way relaying scenario 
by Koike- Akino et al. |2|, and for three-way relaying scenario in 
l3]| . It is observed that adaptively changing the network coding 
map used at the relay according to the channel conditions greatly 
reduces the impact of multiple access interference which occurs 
at the relay during the MA phase. These network coding maps 
are so chosen so that they satisfy a requirement called exclusive 
law. We show that when the four users transmit points from the 
same M-PSK constellation, every such network coding map that 
satisfies the exclusive law can be represented by a 4-fold Latin 
Hyper-Cube of side M. The network code map used by the relay 
for the BC phase is explicitly obtained and is aimed at reducing 
the effect of interference at the MA stage. 

I. Background And Preliminaries 

Physical layer network coding for the two-way relay channel 
exploits the multiple access interference occurring at the relay 
so that the communication between the end nodes can be 
done using a two stage protocol. This two-stage protocol was 
first introduced in f41, and |5|, 16J deal with the information 
theoretic studies for the scheme. In ||2l, modulation schemes to 
be used at the nodes for uncoded transmission for the two-way 
relaying were studied. 

Work done for the relay channels with three or more user 
nodes is given in [31, ifTl- llTol . In Q, a two stage operation 
called joint network and superposition coding, which employs 
four channel uses, three for the MA phase and one for the 
BC phase, has been proposed. The three users transmit to the 
relay node one-by-one in the first phase, and the relay node 
makes two superimposed XOR-ed packets and transmits back 
to the users in the BC phase. The packet from the node with 
the worst channel gain is XOR-ed with the other two packets. 
It is claimed by the authors that this scheme can be extended 
to more than three users as well. Work by Pischella and Ruyet 
in |81 proposes a lattice-based coding scheme combined with 
power control, composed of alternate MA and BC phases, 
consisting of four channel uses for three-way relaying. The 
relay receives an integer linear combination of the symbols 




Fig. 1. A four-way relay channel 

transmitted by the user nodes. It is stated that the scheme 
can be extended to more number of users. These two works 
however deal with the information theoretic aspects of multi- 
way relaying. 

In ifTOl . the authors Jeon et al. adopt an 'opportunistic 
scheduling technique' for physical network coding where us- 
ing a channel norm criterion and a minimum distance criterion, 
users in the MA as well as the BC phase are selected on 
the basis of instantaneous SNR. Their approach utilizes six 
channel uses in case of three-way relaying and it is mentioned 
that the approach can be extended to more number of users. A 
'Latin square-like condition' for the three-way relay channel 
network code is proposed in [9|, and it is suggested that cell 
swapping techniques on these Latin Cubes can be employed to 
improve upon these network codes. The protocol employs five 
channel uses, and the network coding map discussed doesn't 
deal with the channel gains associated with the channels 
explicitly. Latin Cubes have been further explored as a tool to 
find the network coding map used by the relay, depending on 
the channel gain in |3|. The throughput performance of the two 
stage protocol for three-way relaying given in |3 1 is better than 
the throughput performance of the 'opportunistic scheduling 
technique' given in ifTOll at high SNR, as can be observed from 
the plots given in (3]. Our work in this paper further extends 
the approach used in [3J to four-way relaying and employs 
two channel uses for the entire information exchange amongst 
the four users, which makes the throughput performance of 
the scheme better than the other existing schemes 

We consider the four-way wireless relaying scenario shown 
in Fig. 1, where four-way data transfer takes place among the 
nodes A, B, C and D with the help of the relay R assuming 
that the four nodes operate in half-duplex mode. The relaying 
protocol consists of two phases, multiple access (MA) phase, 
consisting of one channel use during which A, B, C and D 




F2XF2= {00.01, 10, 11} 

z<={o, 1,2,3} 
s= { ij,-1,-j} 

Fig. 2. 4-PSK constellation 

transmit to R; and broadcast (BC) phase, in which R transmits 
to A, B, C and D in a single channel use. Network Coding is 
employed at R in such a way that A(/B/C/D) can decode B's, 
C's and D's(/A's, C's and D's /A's, B's and C's) messages, 
given that A(/B/C/D) knows its own message. 

Our physical layer network coding strategy uses a mathe- 
matical structure called a Latin Hyper-Cube defined as follows: 

Definition 1: An n-fold Latin Hyper-Cube L of r-th order 
of side M Ull is an M x M x ... x M (n times) array 
containing M" entries, M""'' of each of M*" kinds, such that 
each symbol occurs at most once for each value taken by each 
dimension of the hyper-cube. 

For our purposes, we use only 4-fold Latin Hyper-Cubes of 
side M on the symbols from the set Zt = {0, 1, 2, t — 1}, 
t > M^. In this hyper-cube, fixing the first dimension, that 
represents A's transmitted symbol, we get a three dimensional 
array, which we will be referring to as cube; then B(/C/D)'s 
transmitted symbols are along the files(/rows/columns) of each 
cube. 

II. Signal Model 

Multiple Access (MA) Phase: 

Suppose A(/B/C/D) wants to send a 2-bit binary tuple to 
B, C and D(/A, C and D/A, B and D/A, B and C). The 
symmetric 4-PSK constellation {±1, ±i} denoted by S as 
shown in Fig. 2, is used at A, B, C and D, and /i : F2 — 5- 5 
denotes the map from bits to complex symbols used at A, 
B, C and D where F2 = {0,1}- Let xa = ^{sa)tXb ~ 
^ (sb) ,xc=tJ. (sc) ,xd — tJ' {sd) G S denote the complex 
symbols transmitted by A, B, C and D respectively, where 
sa,sb,sc,sd G ^2- We assume that the Channel State 
Information (CSI) is not available at the transmitting nodes and 
perfect CSI is available at the receiving nodes. The received 
signal at R in the MA phase is given by. 



Yr = Haxa + Hbxb + Hcxc + Hdxd + Zj, 



(1) 



where Ha, Hb, He and Hd are the fading coefficients 
associated with the A-R, B-R, C-R and D-R link respectively. 
The additive noise Zr ~ CM {Q^a"^), where C7V (O, cr^) 
denotes the circularly symmetric complex Gaussian random 
variable with variance . 

'The definition has been modified slightly from the referred article "On 
Latin and Hyper-Graeco-Latin Cubes and Hyper Cubes" by K. Kishen 
(Current Science, Vol. 11, pp. 98-99, 1942), in accordance with the context. 



Let the effective constellation seen at the relay during the 
MA phase channel use be denoted by Sr {Ha, Hb,Hc, Hb), 
i.e., 

Sii{Ha,Hb,Hc,Hd) = {Haxa + HgXB + Hcxc + HdxdI 

XA,XB,xc,xr) e S}. 

The minimum distance between the points in 
the constellation Sr{HatHbtHc,Hd) denoted by 
dmin{HA,HB,Hc,HD) is given in (O on the next 
page, where 5" = iS x x .. x 5 (n times). From (O, it 
is clear that there exists values of {Ha, Hb, He, Hb), 
for which dmin {Ha, Hb, He, Hd) = 0. Let H ^ 
{{Ha, Hb, He, Hb) G C^ldmin {Ha, Hb, He, Hb) = O}. 
The elements of T-L are called singular fade states. For singular 
fade states, \Sr {Ha, Hb, He, Hd)\ < 4^. 

Definition 2: A fade state {Ha, HB,He, Hb) is defined to 
be a singular fade state for the MA phase of four-way relaying, 
if the cardinality of the signal set Sr {HA,HB,He,HB) is 
less than A"^. 

Let the Maximum Likelihood (ML) estimate of 
{xA,XB,xe,XB) be denoted by {xA,XB,xe,XB) G 5* at R 
based on the received complex number Yr, i.e., 

{xA,XB,xc,xr>) = a.rg min 

(xA ,XB ,xc ,^1 

where H = [Ha Hb He Hd] and X 



\\Yn-HX\\ 



(2) 



[xA xb xe xb] 



Broadcast (BC) Phase: 

During the BC phase, the received signals at A, B, C and 
D are respectively given by, 

111) 

where Xr = j^HaMbMcMd [[xA,XB,xe,XD)) G S' de- 
notes the complex number transmitted by R and H^ , Hg , He , 
and Hjj respectively are the fading coefficients corresponding 
to the links R-A, R-B, R-C and R-D. The additive noises 
Za, Zb, Ze and Zb are CM {{),a'^). During BC phase, R 
transmits a point from a signal set S given by a many to one 
map _\4Ha.Hb.HcMd . 54 _^ ^' chosen by R, depending 
on the values of Ha, Hb, He and Hb- The cardinality of 
S > 64 = 2^, since 6 bits about the other three users needs 
to be conveyed to each of A, B, C and D. 

A cluster is the set of elements in 5^ which are mapped to 
the same signal point in S by the map J^HaMbMcMo ^ Lgj- 
qHaMbMcMo ^ {Ci,C2,..,£.i} denote the set of all such 
clusters. 

Definition 3: The cluster distance between clusters 
Ci, Cj G c"a,Hb,Hc,Hd ^ as given in Q on the next page, 
is the minimum among all the distances calculated between 
the points {xA,XB,xe,XB) G Ci and {x'A,x'B,xe,x'B) G Cj 
in the effective constellation seen at the relay node R. The 
minimum among all the cluster distances among all pairs of 
the clustering C^aMb,Hc,Hd jg jj-g minimum cluster distance, 
as given in (|5]l on the next page. 

The performance during the MA phase depends on the 
minimum cluster distance, while the performance during the 
BC phase is dependent on the minimum distance of the signal 
set iS .A phenomenon referred to as distance shortening, given 
in [4], is described as the significant reduction in the value of 



<^'^in' {Ha,Hb,Hc,Hd) = min \Ha [xa - x'^^) + Hg (xs - x'g) + He {xq - x'^.) + Ho {xd - x'o) | (4) 



dminiHA, Hb,Hc, Hd) = min \Ha (xa - a:'^) + Hb [xb ~ x'g^ + He [xq - x^) + Ho (xd - i-'o) | (3) 

(aiA.^iB ,a:c ,2:0) ,a:J^ )SS* 

(x^ ,X jj')^Cj 

d„,i„ (c"aMb,Hc.Ho'^ = — I Ha {xa - x'^) + Hb {xb - x'g) + He {xc - x'^) + Hd {xd - x'^) | 

M"A.Hb<Hc-Hd (xA,XB,xc,^o)i^M"A.HB.Hc.HD (x\,x'b ,x'^ ,x'^) 

(5) 

rfmin {c'^'^"'^'"'''"'^'"°^\hA,hB,hc,ho) = ^^^i^^^ I hA {xA -x'^)+hB [xb -x'b)+ he {xc - x'^.) + ho {xd " ^d) | 

{x A,x B ,x c ,x o) •(^'a'^'b -^'c '^'d^^'^''' ' 
M^A.Hb.Hc.Hd {xa.xb.xC-^d)T^M"A.Hb.Hc.Hd (x'^,x'b,x'^.x'b,) 
(6) 

j^Ha,HbMc,^d {xa,xb,xc,xd) M^^-^'^-^C'^d [xA,x'g,x'c,x'jy) ,VxA,XB,x'g,xc,x'c,XD,x'j^ eS,{xB,xc,XD) ^ {x'g,x'c,x'j,) 

(V) 

j^HaMb,Hc,Hu {xa,Xb,Xc,Xd) M"^-"''-"<=--"'^ [x'^,XB,x'c,x'^y) ,'iXA,x'A,XB,Xc,x'c,XD,x'g, eS,{xA,Xc,XD) + {x'j^,x'c,x'o) 

(8) 

j^Ha,Hb,Hc,Ho (xA,XB,xc,xu)i= M"^'"b'"^-"° {x'j^,x'g,xc,x'j^) ,VxA,x'^,XB,x'g,xc,XD,x'^j a S,{xa,xb,xd) + {x'^,x'g,x'j^) 

(9) 

j^Ha,Hb.Hc,Hd {xa,xb,xc,xd) + m"^-"''-"c-Hd (x^, x'^, z^, xo) , Vxa, x^^, zs, x^, xc, x'^ £ S,{xa,xb,xc) ^ (x'j^,x'g,x'(.) 

(10) 



d„^^n [c" A,H b Mc .H for yalues of {H A. H B , He , H d) in 
the neighborhood of the singular fade states. If the clustering 
used at the relay node R in the BC phase is chosen such 
that the minimum cluster distance at the singular fade state is 
non zero and is also maximized, then the effect of distance 
shortening can be avoided. 

We say a clustering C^a,Hb,Hc.Hd removes a singular fade 
state {Ha,Hb,Hc,Hd) G H, if d„„„ (C^^-^«^^^'^°) > 
0. Let c^{Ha,HbMcMd)} denote the clustering which re- 
moves the singular fade state {Ha, Hb, He, Hb) (selecting 
one randomly if there are multiple clusterings which remove 
the same singular fade state {Ha, Hb, He, Hb))- Let the 
set of all such clusterings be denoted by C-u, i.e., C-h = 
{Ci^HA.HB.HcMn)} . {H A, H b , He , H b) e n). 

Definition 4: The minimum cluster distance of the cluster- 
ing c'^(Ha,Hb.Hc,Ho)} foj {HA,HB,He,HB) e H, when 
the fade state {Ha, hs, he, ho) occurs in the MA phase, 
denoted by d„„„ {c{{Ha,Hb,Hc,Ho)} ^ he, he, ho), is the 
minimum among all its cluster distances. 

When {hA,hB,hc,hB) ^ 'H, we choose the 

clustering C^A,hB,hc,hn ^ C^iHAMB,Hc.Hn)} c-H, 

such that d„,,^{C^("A,HBMc,Hn)}^hA,hB,hc,hD) > 

drmn {cH^A'^B '"'c 'H'^)} , Ha, h g , he , ho) ,y (H A, H g , He , H o) + 

{H'j^,H'g,H'^,H'j^ G H. The clustering used by the relay is 
indicated to A, B, C and D using overhead bits. 

In order to ensure that A(/B/C/D) is able to decode B's, C's 
and D's(/A's, C's and D's /A's, B's and D's /As, B's and C's) 
message, the clustering C should satisfy the exclusive law, 
as given in (|7|l, ©, ©, (fTOl i on the next page. We explain 
Exclusive Law in more detail in the next section. 
The contributions of this paper are as follows: 
• Exchange of information in the wireless four-way relay- 
ing scenario is made possible with totally two channel 



uses using our proposed scheme. 

> It is shown that if the four users A, B, C, D transmit points 
from the same M-PSK constellation, the requirement of 
satisfying the exclusive law is same as the clustering 
being represented by a 4-fold Latin Hyper-Cube of side 
M. (Section III) 

* The singular fade states for the four-way relaying scenario 
are identified. (Section IV) 

> Clusterings that remove these singular fade states are 
obtained and are of size between 64 to 90. (Section V) 

« Simulation results are provided to verify that the adaptive 
clustering as obtained in the paper indeed performs better 
than non-adaptive clustering at high SNR. (Section VI) 

The remaining content is organized as follows: Section 
III demonstrates how a 4-fold Latin Hyper-Cube of side 4 
can be utilized to represent the network code for four user 
communication. In Section IV we define what we mean by 
singular fade subspaces and in Section V, focus in on the 
removal of such singular fade subspaces using 4-fold Latin 
Hyper-Cube of side 4. Simulation results are provided in 
Section VI. Section VII concludes the paper 

III. The Exclusive Law and Latin Hyper-Cubes 

In the previous section, it was stated that in order to ensure 
that A(/B/C/D) is able to decode B's, C's and D's(/A's, C's and 
D's /A's, B's and D's /A's, B's and C's) message, the clustering 
C that represents the map used at the relay should satisfy 
the exclusive law. We assume that the nodes A, B, C and 
D transmit symbols from the 4-PSK constellation. Consider 
a4x4x4x4 array, containing 4^ — 256 entries indexed 
by {xa,xb,xc,^d), i-S-, the four symbols sent by A, B, C 
and D in the MA phase. The four cubes of this 4x4x4x4 
array, are indexed by four values taken by xa- Each file (/row 




Fig. 3. A 4-fold Latin Hyper-Cube of side 4 represents the exclusive law 
constraint for the relay map when 4-PSK is used at end nodes 

/column) of each such cube is indexed by a value of xb {Ixc 
Ixd), for this fixed value of xa- The repetition of a symbol 
in a cube results in the failure of exclusive law given by 
Considering the 4x4x4 array with its files being the 
first(/second/third/forth) files of the 4 x 4 x 4 x 4 aiTay. Each 
4x4x4 array so obtained, corresponds to a single value of 
Xb- a repetition of a symbol in this array will result in the 
failure of exclusive law given by ([8j. Similarly, a repetition of 
a symbol in the 4x4x4 array with its rows(columns) being 
the first(/second/third/forth) columns of the 4x4x4x4 array, 
that corresponds to a single value of xc{xd), will result in 
the failure of exclusive law given by (l9]l((fT0]i). Thus, for the 
exclusive law to be satisfied, the cells of this array should be 
filled such that the 4x4x4x4 array so obtained, is a 4- 
fold Latin Hyper-Cube of side 4, for i > 64 (Definition 1). 
The clusters are obtained by putting together all the tuples 
{i,j,k,l),i,j,k,l e 0,1,. ..t — 1 such that the entry in the 
k, l)-th slot is the same entry from Zt. 
The same can be extended to an M-PSK constellation, 
i.e., if the network code used by the relay node in the 
BC phase satisfies exclusive law, it can be represented by 
a 4-fold Latin Hyper-Cube of side M and t > with 
cubes(files/rows/columns) being indexed by the constellation 
point used by A(/B/C/D), symbols from the set Zm (Fig. 3). 
Any arbitrary but unique symbol from {Ci, ■■,Ct} denotes a 
unique cluster of a particular clustering. 

IV. Singular Fade Subspaces 
As stated in Section II, a clustering q^a^Hb-HcHd jg g^id 
to remove singular fade state {Ha, Hb, He, Hu) e if 
dmin (C^a-^bMcMd"^ > [ ^jjy j-yyo mcssagc sequences 
{xA, xb,xc, xd) £ that coincide in the effective constella- 
tion received at the relay during the MA phase should be kept 
in the same cluster of C^a,Hb,Hc,Hd ^ removing singular 
fade states for a four-way relay channel can alternatively be 
defined as: 

Definition 5: A clustering c"a,HbMcMd jg ^^^^ ^g. 
wove the singular fade state (Ha, Hb, He, Ho) G T-L, if 
any two possibilities of the messages sent by the users 
{xa,xb,xc,xd) ,{x'j^,x'g,x'u, x'jj ) G iS^ that satisfy 

Haxa + Hbxb+Hcxc + Hbxd = H ax'a+ H Bx'g+ Hc.x'c + H Dx'jj 

are placed together in the same cluster by the clustering. 

Definition 6: A set {{xatXbtXCtXd)} G S'^ consisting of 
aU the possibilities of {xa,xb,xc,xo) that must be placed 
in the same cluster of the clustering used at relay node R 



in the BC phase in order to remove the singular fade state 
{Ha, Hb, He, H]j) is referred to as a Singularity Removal 
Constraint for the fade state {Ha, Hb,Hc, Hd) for four-way 
relaying scenario. 

The relay receives a complex number, given by ([T]i, at the 
end of MA phase. Using the ML estimate of this received 
complex number, R transmits a point from the constellation 
S' with cardinality at most A"^. Instead of R transmitting a 
point from the 4** point constellation resulting from all the 
possibilities of {xatXb,xc,xd), depending on the fade states, 
the relay R can choose to group these possibilities into clusters 
represented by a smaller constellation, so that the minimum 
cluster distance is minimized, as well as all the users receive 
the messages from the other three users, i.e., the clustering 
satisfies the exclusive law. We provide one such clustering for 
the case of four-way relaying in the following. 

Let {Ha, Hb, He, H]j) be the fade coefficient in the 
MA phase. Suppose {Ha, Hb, He, Hu) is a singular fade 
state, and F is a singularity removal constraint corre- 
sponding to the singular fade state {Ha, Hb, He, Hd). 
Then there exist {xA,XB,xc,XD),{x'j^,x'g,x'fj,x'jj) £ F, 
{xa,xb,xc,xd) ^ {x'j^,x'g,x'(j,x'jy) such that. 



+ Hb^b 
■'a) + Hi 



{XB> - x' ) = 



>(Ha.Hb.Hc.Hb,) e 



where for a 4 x 1 non-zero vector v over C, 

(v)'^ — < 10 — {wi, 102, 103 , W4) E \ WiVi -\- W2V2 + W3V3 + W^V^ ~ ^ I 

j_ (13) 

It can be proven that {v) is a three-dimensional vector 
subspace of C^. Since XA,XB,xc,XD,x'j^,x'g,XQ,x'jj € S, 
where S is finite, there are only finitely many possibilities 
for the right-hand side of (fT2] i. Thus the uncountably infinite 
singular fade states {Ha, Hb, He, Hb), are points in a finite 
number of vector subspaces of C*. We shall refer to these finite 
number of vector subspaces as the Singular Fade Subspaces 
[3 1. There are four possibilities of singular fade subspaces for 
four-way relaying as we explain individually in the following. 
Case 1: Only one of xa~ x'j^, xb ~ x'g, xc — x'q or xd — x'jy 
is non-zero, resulting in the following 4 subcases: 

1) Xb — x'g, Xc = x'fj, Xd = x'j^ and xa 7^ x'j^ 

2) XA = x'^, Xc 

3) XA — x\, Xb 

4) XA = x'j^, XB 
Case 2: Two of xa - 

non-zero, resulting in the following 6 subcases: 



= x'(-,, xb) = x']j and xb ^ x'g 

- x'g, Xb — x']j and xc x'^ 

- x'g, Xc = x'q and xb ^ x'jj 
x'a^ Xb — x'g, Xc — x'(j or Xb — X 



are 



1) 


XA 


— ^Ai 


Xb 


— ^'b. 


Xc ^ x'c 


and Xd 7^ x'j~, 


2) 


XA 


— ^A-> 


Xc 


— 


XB ^ x'g 


and Xd x'^ 


3) 


XA 


= Xj^, 


Xd 




XB 7^ x'g 


and Xc 7= x'^ 


4) 


Xb 


— Xb, 


Xc 


= Xc, 


XA 7^ X 


and Xd x'^ 


5) 


Xb 


— ^B, 


Xd 


— Xd, 


XA 7^ x'^ 


and Xc ^ a^'c 


6) 


Xc 


— 


Xd 


— ^D, 


XA 7^ x'j^ 


and Xb ^ x'q 



Case 3: Three of xa — x'j^, xb — x'g, xc — x'^-, or xd — x 
are non-zero, resulting in the following 4 subcases: 



D 



1) XA = x'j^, XB 7^ x'g, xc ^ x'(j and xd ^ x'j^ 

2) XA ^ x'j^, XB = x'g, Xc 7^ x'(j and xd ^ x'^ 

3) 7^ ^'a' 7^ x'g, Xc = x'(j and xd a;^ 

4) ^ x^, ^ x^, Xc ^ x'(j and x_d = x^ 

Cflse 4; All of xa ~ x'a, xb — x'g, xc — x'q or xb — x^ are 
non-zero, i.e., xa x^, xs 7^ x^, xc 7^ Xp and xd 7^ x'^. 

Case i; Only one of xa — x^, xb — x^, xc — x'q ot xb ~ x^ 
is non-zero, resulting in 4 subcases. Without loss of generality, 
we discuss the first subcase of Case 1. The singular fade sub- 



space in this case will be given by S' 



XA ■ 



Let V denote the set of differences of the points of 5, i.e., 
V={xi- Xj I Xi, Xj G 5} = {0, ±1 ± J, ±2j, ±2} . We can 
write, V = {0} U I?i U V2, where X»i = {±1 ± j} and 
V2 - {±2j, ±2}. 



Now, XA - 



G V can take eight non-zero values. So there 



are eight total possibilities for the vector [xa — x^ , 0, 0, 0] , 
where v* denotes the transpose of a vector v. Since each one 
of ±1 ± j, ±2j and ±2 can be obtained as scalar multiples 
of 1 + j (over C), we have. 



0, 0, 0]^: 

0, 0, 0]^: 

+ 0, 0, 0]^ : 

+ 0, 0, 0]^ : 

0, 0, 0]^: 

0, 0, 0]^: 

0, 0, 0]^ : 



l-l-j, 0, 0, 0]^ 
0, 0, 0]^ 

0, 0, of 

I2j, 0, 0, of 
[-2i, 0, 0, 0]^ 
[2, 0, 0, 0]^ 
[-2, 0, 0, 0]^ 



Therefore, 



" 1 + J " 




r ±1 ± J 1 




r 1 




' ±2 ■ 





)< 





)< 















n 









n 







n 









. n - 



As a result, only one singular fade subspace arises from 
this subcase, viz., 5' = <^[l + j, 0, 0, 0]^)^ . Similarly, for the 
other three subcases, there is a single singular fade subspace, 
resulting in a total of 4 singular fade subspaces arising from 
this case. 



Xb — x'g, Xc — x'q or Xd ~ x^ are 



for this case is given by S" 



Both 



Case 2: Two of xa — x 
non-zero, resulting in 6 subcases. Without loss of generality, 
we consider the case for which the singular fade subspace 

Xa — x'j^ 

Xb - x'g 



XA — x'j^ and xb — ^'b ^ ^ '^^'^ '^^^e eight non-zero values 
each. There are therefore, 64 total possibilities for the vector 

[xa ~ X'j^, XB - X'g, 0, 0]*. 

Lemma 1: Given v = \^xa — a;^, xb — x'^, xc — x'^, xo — a^'o] 
with entries from 2?i U V2, when at least one of xa ^ 
3^A' "^B ^ x'g, Xc 7^ x'(j and xd / x^, there are precisely 
4 or 8 vectors (including v) with entries from Vi U that 
generate the same vector space over C as v. 



We partition and prove this Lemma in three parts: Lemma 
2, Lemma 3, Lemma 4 given in Case 2, Case 3, and Case 4 
respectively of this paper. 

Lemma 2: When xa 7^ x^, xb 7^ x^, xc — 
x'(j and Xd = x'jj, for a given vector v — 
[xa — x'j^, Xb — x'g, 0, 0]* over Vi U 2?2, there are 
precisely 4 or 8 vectors (including v) over T>i U T>2 that 
generate the same vector space over C as v. 

Proof: The difference constellation V — AS — 
{s — s' : s, s' G S} where S is 4-PSK constellation, is given 
by m, 

AS = {0} U {2sm (7m/4) e^'^^/^ln odd} 

U {2sm(7rn/4)e^(''''/2+V4)|^ g^gnj ^ 

where 1 < n < 2 and < fc < 3. So, v can also be written 
as, 

2sin^e^'t'^ 

2sin^e^'l'^ 






XA 


— x'^ 




Xb 


~ Xg 



















where 
even. 



kiTT/2 if ki is odd and (f>i — kiTr/2 + 7r/4 if fc; is 



A vector w over 2?i U shall generate the same vector 
space over C iff uj is a scalar multiple of v, i.e., for some 
complex number reJ^ G C, 

2sin^e^'^^ ' 

2sin^ei'l>^ 




where for i — 3,4 (/>i = 
(f)i = kiTT/2 + tt/A if hi is even. 



2sin^e^'f'^ 
2sin^e''*'* 






fci7r/2 if ki is odd and 



Then, 



4 4 



and 



7rfc2 

2sin e 

4 



092 



4 



(14) 



(15) 



Dividing ( fT4] i by (fTsT i and taking modulus of both sides, we 



get 



sm 



Trki 



Sin 



(16) 

4 ■J""' 4 ■ 

As shown in 1121 . this is possible only if ki = and k2 — k^. 
Also, from ( fT4b and (fTSl) we have 



.7. ^ 



sm- 



.7., 



1) 



From ( fTSI l and ( fTTI ). we have 



(17) 



(18) 



Fig. 4. Null Spaces of the Singular Fades Subspaces for 

Note that here, the LHS is fixed. It therefore suffices to 
compute the number of values that RHS takes for the fixed 
value of the LHS. It can be verified, that for a fixed value of 
(j^i — 4'2, there are precisely four pair of values of (p^ and (/)4 
that result in the same value of 03 — 04. We now look at the 
following possibilities: 

Case 1: ki = k2- Then, — k^, i.e., there are exactly two 
possibilities for fci and k2, viz., /ci = ^2 = 1 and ki = k-z = 2. 
With two pairs of values for k^ and /c4 and four pairs of values 
for 03 and 04, we have a total of eight set of values that w 
can take. Hence, in this case, the vector space generated by v 
can be generated by exactly eight other vectors over Pi U2?2- 
Case 2: ki ^ k2- Then, k^ ^ k^, i.e., there is precisely one 
possibility for fci and fc2, viz., fci = k^ and ki — k^. With 
only one possible set of values for k^ and kj^ and four pairs of 
values for 03 and 04, we have a total of four set of values that 
w can take. Hence, in this case, the vector space generated by 
V can be generated by exactly four other vectors over I?i Ul?2- 

■ 

In this case we end up with 12 singular fade subspaces 
given by the null spaces of the spaces given above in Fig. 
m Similarly, for each of the other five subcases, there are 12 
singular fade subspaces, resulting in a total of 72 singular fade 
subspaces for the case. 



the case ^a' ^b' 



and X£i = X , 



Case 3: Three of xa 



XB 



XC 



C q or xd 



are non-zero, resulting in 4 subcases. Consider the subcase, 
where the singular fade subspace in this case is given by, 

xa- x'j^ ^ 



x'g and Xc 



Xg - X' 

xc - x[ 


— Xq 



generate the same vector space over C as v. 

Proof: The proof of this Lemma is similar to the proof 
of Lemma 1 and is therefore omitted. ■ 

Since xa — x'^^, xb — x'g and xc — Xq ^ V are non-zero, 
we can say that 

Xa — x'j^, Xb — Xg and xc — x'c G Pi U P2- 

As a result, we have the following three subcases: 

1) One of xa — a;^, xb — x'^ and xc — x'c € Pi 



2) Two of XA 



^A' 



Xb — x'g and xc — x'^j G Pi 



3) All of Xa — x'j^, Xb — x'g and xc — x'^j e Pi 
We deal with each one of the subcases one-by-one. 
Subcase 1: One of xa — x'^^, xb — x'g and xc — x'c 
Without loss of generality, we assume that xa — x'j^ 



e Pi. 

e Pi 



and Xb 



Xc 



e P2. There are 64 possibilities 



for the vector v' = [xa^x'j^, xb—x'^, xc—x'q, 0]*. 
For each of the 64 possibilities for the vector v' — 
[xA—x'j^, Xb ~ x'g, Xc — x'c, 0]*, precisely 4 vectors of 
length 4 over Pi U P2 generate the same vector space over 
C. As a result, this case leads to 16 singular fade subspaces. 
The same holds for the case when xb — x'g e Pi and 



XA 
XA 



Xc 

Xb - 



£ P2, or when xc 



G Pi and 



S P2- This subcase therefore results in 



48 singular fade subspaces. 
Subcase 2: Two of xa 



Xb — x'g and xc 



e Pi. 



This subcase can be dealt with similar to the Subcase 1 and 
hence also results in 48 singular fade subspaces. 



Subcase 3: All of xa 



Here each of xa 



Xb — x'g and xc 



X 



e P can take eight non-zero values. 



A, Xb — Pi (or P2). There are 128 possibilities for the vector 

[xa — x'j^, Xb — x'g, Xc — x'^, 0]* over Pi (or P2). For a 



There are therefore, 512 total possibilities for the vector 

[xA-x'j^^, XB-x'g, xc-x'(j, 0]*. 

Lemma 3: For the case when xa ~ x'^^ 0, xb 
0, Xc — x'c and xd — x'r 



X 



: 0, for a given vector 
Xb — x'g, Xc ~ x'(j, 0]* over Pi UP2, there 
are precisely 4 or 8 vectors (including v) over Pi U P2 that 



[XA 



given such vector, 8 other possibilities of 4 length vectors 
over Pi U P2 that generate the same vector space over C. The 
subcase therefore leads to a total of 128/8=16 singular fade 
subspaces. 

The three subcases result in a total of 48-1-48-1-16=112 
singular fade subspaces. So, there are a total of 4 x 112 = 448 
singular fade subspaces corresponding to Case 3. 



Case 4: All of xa — a;^, xb ^ x'g, xc — Xq or xd — x'q 
are non-zero. The singular fade subspace in this case is given 

/ r ^A-x'^ 1 \ ^ 

by, 5""=/ x-c-^^ / • Here each of -a:^^, - 

\ L ^d-x'd J / 
x'bj xc — x'(j and xd — x'j^ G V can take eight non-zero 

values. There are therefore, 4096 total possibilities for the 

vector [xA — x'a-, xb — x'g, xc — x'fj, xo — x'^]* . 

Lemma 4: For the case when xa ^ x'j^, xb 7^ 

x'bi Xc 7^ x'^ and xo ^ x^, for a given vector v ~ 

[xA - x'j^, xb - x'g, Xc - x'fj, xd - x'j^f over Vi U V2, 

there are precisely 4 or 8 vectors (including v) over 2?i U 1^2 

that generate the same vector space over C as v. 

Proof: As mentioned in the proof of Lemma 1, from lfT2l . 

V = AS = {0}U [2sin (7m/4) e^''''/^\n odd} 

U {2sm(7rn/4) e^'(''''^/2+7r/4) |^ g^gj^j ^ 

where 1 < n < 2 and < fc < 3. Therefore, we can write, 

4 ^ 









2sin 




— X^ 








— X^ 




2sin 


Xc 


— X(j 




2sin 


X_D 


— X]j 










2sin 



where 
even. 



ki'K/2 if ki is odd and 0,; — ki'K/2 + 7r/4 if fcj is 



A vector w over Vi U ^2 shall generate the same vector 
space over C iff w is a scalar multiple of v, i.e. for some 
complex number re^^ G C, 



2sin^e^*^ 
2sin^ei'f'-'- 



2sin^^e 



4 



2sin^e't'^ 



where for i = 4, 5,6 1/)^ = ki-K/2 if fc^ is odd and 
(pi = ki-K/2 + 7r/4 if fc^ is even. 



Then, 



2sin e^^^ — re-' x 2sm e-' 

4 4 



o . 7rfc2 

isin e 

4 

2sm — ^6' 
4 



4 



_ ^gj-e X 2sin e 

4 



J 07 



and. 



Zsin e-' 

4 



re^^ X 2sin—e^'f'' 
4 



(19) 
(20) 
(21) 

(22) 



Dividing ( fT9] l by (|20] | and taking modulus of both sides, we 



get 



4 



(23) 



As shown in 1 12], this is possible only if fci = k^ and k2 — k^. 
Similarly, dividing ( fT9] l by ( |2TI ) and taking modulus of both 
sides, we get 



sm 



sm 



(24) 



This is possible only if fci = ^5 and k^ = k^. Similarly, 

ki = fcg and ki = ks- (25) 
Also, from ( fT9l ) and ( |20l ) we have 



sm 



TTke 



4 4 

From ( [23T l and ( [26b . we have 



Similarly, 



and 



gj(01-02) _ gj(05-06)_ 



fe) — gj(05-07) 



1) 



(26) 



(27) 



(28) 



(29) 



In (Elll, (ESll and (l29]l, the LHS is fixed. It therefore suffices 
to compute the number of values that the RHS takes for fixed 
LHS in the three equations. It can be verified, that for a fixed 
value of (pi, (j)2, (f>3 and (j)4, there are precisely four set of 
values of 05, ^g, 07 and (ps that result in the same value of 
01 — 02, 01 — 03 and 01 — 04. We now look at the following 
possibilities: 

Case 1: fci = ^2 = ^3 = ^4. Then, k^ ~ k^ ~ ky = fcg, i.e., 
there are exactly two possibilities for fcs, fcg, fcy and fcg, viz., 

fcs = fce = fcy = fcs = 1 or fcs = fcg = fc? = fcs = 2. With two 
sets of values for fc5, k^, ky and fcg and four sets of values for 
05, (pQ, 07 and 08, we have a total of eight set of values that w 
can take. Hence, in this case, the vector space generated by v 
can be generated by exactly eight other vectors over Pi UX'2- 
Case 2: At least one of fci 7^ for 1 < i < 4, 1 < j < 4. 



Then, fci+4 7^ fcj+4, so that, using (1231 1. (l24] | and (l25l l. there 
is precisely one possibility for fcs, fcg, fc7, fcs, viz., k^ = ki, 
fce — ^2, fc? = fcs and fcg = fc4. With only one possible set 
of values for k^, fcg, fcy and fcg and four sets of values for 
05, 0g, 07 and 08, we have a total of four set of values that w 
can take. Hence, in this case, the vector space generated by v 
can be generated by exactly four other vectors over Vi U 2?2- 

■ 

In Case 4, since all of XA—a;A I xb—x'^, xc~x'c and xd — 
x'jj 6 T> are non-zero, we can say that 



Xa ~ X A^ Xb ~ xq, Xc 



and Xd — x'^, G 2?i U 2?2 ■ 



As a result, we have the following four subcases: 



1 ) Only one of xa 
Vi 

2) Two of Xa 



c^, XB—x'g, xc'—x'fj and x^ — x'^ G 



3) Three of Xa— a;' 

'X^, Xs 



4) AllofxA- 



xb- 

A> 



'X^, xc~X(j 

-Xg, XC~Xq 

Xb^ xc~-^c 



and xd — x'jj G Pi 
and X£)— x'r, G Pi 



xj-, and X£) — x'r. 



ePi 



We deal with each one of the subcases one-by-one. 



Subcase 1: Only one of xa — x'^^, xb — x'g, xc — 
x'q and xb — x'j~, € 2?i. Without loss of generality, we 
assume that xa — x'^ £ Vi and xb — x'g, xc — x'q, xd — 
x'jj G 1^2. There are 256 possibilities for the vector v' = 
[xa — x\, Xb — x'q, Xc — x'q, xd — x'j^f . But precisely 4 
vectors of length 4 over T>i i-e. the {±1, ±j} scalar mul- 
tiples of the vector generate the same vector space over C. As 
a result, the case xa—x'j;^ € Vi and XB—x'g, xc—x'qG V2 
leads to 64 singular fade subspaces. The same holds for the 
case when XB—x'g £ Vi and XA — x'j^^, xc—x'q, xd—x'jj& 
'D2,whenxc — x'u ^Vi, and XA—x'j^, xb—x'^, xd—x'j^^ 
X>2, or when x/)— £ 2?i x^— x^, xb—x'^, xd—x'jj € 2?2- 
This subcase therefore results in 256 singular fade subspaces. 

Subcase 2: Two of xa—x'^, xb —x'g, xc — x'q and xd — 
x'jj e Vi. This subcase can be dealt with similar to the 
Subcase 1 and results in 384 singular fade subspaces. 

Subcase 3: Three of X A — x'j^, xb—x'^, xc — x'q andxD — 
x'jj &Vi. This subcase can also be dealt with similar to the 
Subcase 1 and results in 256 singular fade subspaces. 

Subcase 4: All of, or none of xa — x'j^, xb — x'g, xc — 
x'q and Xd — x'jj e X>i. There are 256 possibilities for 
the vector [xa — x'^, xb — x'^, xc — x'q]* over Vi. For a 
given such vector, possibilities of other 4 length vectors over 
T>i U T>2 that generate the same vector space over C are 
the {±1, ±j, ±1 ± j'} scalar multiples of the vector. The case 
Xa — a;^, xb — x'g, xc — x'q € Vi leads to a total of 64 
singular fade subspaces. These 64 singular fade subspaces have 
been listed in the Appendix. 

The four cases result in a total of 4+72+448+960=1484 
singular fade subspaces. As done for 4-PSK case, this com- 
putation can be similarly done for the case when M-PSK 
(M > 4) constellation is used at the nodes. We now discuss 
how these singular fade subspaces can be removed using 4- 
fold Latin Hyper-Cube of side 4. 

V. Removing singular fade subspaces 

We now cluster the possibilities of {xa,xb, xcxo) into a 
clustering that can be represented by a 4-fold Latin Hyper- 
Cubes of side 4. This clustering is represented by a con- 
stellation given by S' , which is utilized by the relay node 
R in the BC phase. The objective is to obtain a clustering 
that removes the singular fade subspaces, and then attempts 
to minimize the size of this constellation used by R. This 
clustering to be used at R, first constrains the possibilities of 
[xa, xb,xc,xd) received at the MA phase, with the objective 
of removing the singular fade subspaces, and fills the entries 
ofa4x4x4x4 array representing the map to be used at the 
relay using these constraints, and then completes the partially 
filled array obtained to form a 4-fold Latin hyper-cube of side 
4. The mapping to be used at R can be obtained from the 
complete Latin hyper-cube keeping in mind the equivalence 
between the relay map with the 4-fold Latin Hyper-Cube of 
side 4 as shown in Section III. 

During the MA phase for the four-way relaying scenario, 
nodes A, B, C and D transmit to the relay R. As shown 
in the previous section, there is a total of 1484 singular 



fade subspaces. Let the fade state (/ia, ^b, ^c, ^d) denote 
a point in one of the 1484 singular fade subspaces. The 
constraints on the 4x4x4x4 array representing the map 
at the relay node R during BC phase for a singular fade 
state, can be obtained using the vectors of differences, viz.. 



[XA -X\Xb- X'g, XQ - x'q, Xd 



contributing to 



this particular singular fade state. So, if {h,A,hB,hc,hD) € 
[xa — x'j^ Xb — x'g, xq — x'q, Xd — x'jj], then, for 

{xa,Xb,XC,Xd) , {x'j^,x'q,x'q,x'jj) e S'^, hAXA + hBXB + 

hcxc + hDXD = hAx'jn + hBx'g + Hcx'q + hax'^. 
For a clustering to remove the singular fade state 
{hA,hB,hc,hD), i.e., for the minimum distance of the 
clustering to be greater than (Section 11), the pair 
{xa,xb,xc,xd) ,{x'A,x'g,x'Q,x'jj) must be kept in the 
same cluster. Alternatively, we can say that the entry 
corresponding to {xa,xb,xc,xd) in the 4 x 4 x 4 x 4 
array must be the same as the entry corresponding to 
{x'^,x'g,x'Q,x'j)). Similarly, every other such pair in <S^ 
contributing to this same singular fade subspace must be kept 
in the same cluster. Apart from all such pairs in <S^ being 
kept in the same cluster of the clustering, in order to remove 
this particular fade state, there are no other constraints. This 
constrained 4x4x4x4 array can then be completed by 
simply fining the top-most and the left-most empty cell in 
the earliest file of the earliest cube, with Ci,i > 1 in the 
increasing order of i such that the completed array is a 
4-fold Latin Hyper-Cube of side 4. This algorithm is given in 
detail in Algorithm 1 . The above clustering scheme, however, 
cannot be utiUzed to remove the singular fade subspaces of 
Case 1, Case 2 and Case 3 of the previous section, as shown 
in the following lemma. 

Lemma 5: The clustering map used at the relay node R 
cannot remove a singular fade state corresponding to the 
Case 1, Case 2 and Case 3 of the previous section and 
simultaneously satisfy the mutually exclusive law. 

We refer to the singular fade subspaces given in Case 1, 
Case 2 and Case 3 of the previous section, whose harmful 
effects cannot be removed by a proper choice of the cluster- 
ing, as the non-removable singular fade subspaces. We now 
illustrate the removal of a Case 4 singular fade state with the 
help of the following example. 

Example 1: Consider the case for which the singular fade 
subspace is given by 



-1 + J 
1 + 3 
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The first vector is l+j, 1+i, 1— i]- Now, 

— 1 — j can be obtained either as a difference of xa = — 1 
and x^ = j or as a difference of xa = —j and x'^ = 1. 
Similarly, I -\- j can be obtained as a difference of 1 and 
—j or as a difference of j and —1; 1 — j can be obtained 
as a difference of xd = 1 and x^ = j or as a difference 
of Xd = —j and x^ = 1. Thus, the entries corresponding 
to 1, 1, 1), {j, -j, -j,j)}, {(-1, 1, 1, -j), (j, -j, -j, -1)}, 

{(-1, 1), {j, -j, {(-1, -j), {j, -j, -1, -1)}, 



-1)}. 
-1)}. 
-1)}. 
-1)}. 
-1)}. 

must 



{(-1, j, 1, 1), u, -1, {(-1, j, 1, -j), U, -1, , - 
1), (i, -1, j, -i), (i, -i, - 

{{-j, 1, 1, 1), (1, -j, 1, 1, (1, -j, - 
{(-i, i,j, 1), (1, -i, -1, j)}, {{-j, -j), (1, -i, - 

1, 1), (1, -1, -j,m, 1, -i), (1, -1, - 

1), (1, -1, -1, j)}, -i), (1, -1, -1, -1)} 

be the same in the 4x4x4x4 array representing the 
clustering. Similarly the other constraints can be obtained. 
Replacing —1, — j with 0, 1, 2, 3, these constraints can be 
represented in 4x4x4x4 array representing the clustering. 
The constrained Hyper-Cube of side 4 is completed using 
Algorithm 1 to form a 4-fold Latin Hyper-Cube of side 4. 
The completed Latin Hyper Cube is as shown in Table J] 

Similarly, for each of the 960 possibilities of singular fade 
subspaces of Case 4 of the previous section, a clustering of 
size between 64 to 90 can be achieved by first constraining 
the array representing the relay map in order to remove the 
singular fade state and then completing the constrained array 
using the provided algorithm. 

Algorithm 1: Obtaining the 4-fold Latin Hyper-Cube of 
side 4 from the constrained 4x4x4x4 array 
Input: The constrained 4x4x4x4 array 
Output: A 4-fold Latin Hyper-Cube of side 4 

representing the clustering map at the relay 

1 Start with the constrained 4x4x4x4 aiTay X 

2 Initialize all empty cells of A:" to 

3 The k, element of X is the i*'' cube, the j*'* 
file, the k*^ row and the column cell. 

4 for 1 < i < 4 do 

5 for 1 < j < 4 do 

6 for 1 < fc < 4 do 

7 for 1 < Z < 4 do 
if cell {i, j, fc, /) of X is NULL then 

9 Initialize c=l 

10 if Cc does not occur in the i^^ cube 
of X, the j*'* file of X, the fc*'' row of 
X and the l*-^ column of X then 

11 replace at cell (i, j, fc, I) of X 
with Cc, 

else 

I c=cH-l; 
end 

end 



12 
13 
14 
15 
16 
17 
18 

19 end 



end 



end 



end 



ensuring a minimum cluster distance greater than zero for 
the fade state. This scheme removes a subset of singular fade 
states, called removable singular fade states. Simulation results 
presented in this section, done for the case when the end nodes 
use 4-PSK signal set, identify some cases where the proposed 
scheme outperforms the naive approach that uses the same 
map for all fade states and vice verse. 
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Fig. 5. SNR Vs BER curves for different schemes 



Here, Ha, Hb, He, Hd, H'^, H'g, H'^^ and H'j^ are dis- 
tributed according to Rician distribution and channel variances 
equal to dB. The plot for the case with a Rician Factor of 
20 dB for a frame length of 256 bits is as shown Fig. |5] The 
figure shows the SNR vs bit-error-rate curves for the following 
schemes: (a) the adaptive network coding scheme presented 
in this paper; (b) the non-adaptive network coding using two 
channel uses, in which the same 4x4x4x4 array is used by 
relay as an encoder for all channel conditions, given by Table 
nil It can be seen from Fig. |5] that the scheme based on the 
adaptive clustering relaying perform better than the schemes 
based on non-adaptive clustering at high SNR, since adaptive 
clustering removes 960 singular fade states. 
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Adaptive Network Coding for four-way relaying using two 
dianncl uses when 4-PSK is used at end nodes 
^ Non-Adaptive Network Coding for four-way relaying using two 
channel uses when 4-PSK is used at end nodes 
Non- Adaptive Network Coding for four-way relaying using tliree 








channel uses when 4-P5K is used at end nodes 
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Fig. 6. SNR Vs Throughput curves for different schemes 



VI. SIMULATION RESULTS 

Our scheme enables information exchange in the four way 
relaying scenario amongst the four users in total two channel 
uses. It is based on the removal of singular fade states for 
the case. Removing a singular fade state essentially means 



Fig. |6l shows the throughput comparison between (a), (b), 
and (c): a non-adaptive network coding scheme that uses three 
channel uses, viz.. A, B transmit in the first channel use, C, D 
transmit in the next channel use, and the relay node transmits 
the function of these MA phase transmissions in the third 
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Example 1: Latin-Hyper Cube representing the Relay Map where entries xa's and xb's are as mentioned, xc's entries are along 
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TABLE n 

The 4-fold Latin Hyper-Cube of side 4 used at the relay node as an encoder in the BC for non- adaptive network coding using two 

channel uses for all channel conditions 



channel use. At high SNR regime, the scheme presented in 
this paper leads to higher throughput, i.e., (a) outperforms both 
(b) and (c). 

VII. Conclusion 

We considered the four-way wireless relaying scenario, 
where four nodes operate in half-duplex mode and trans- 
mit points from the same M-PSK constellation. Information 
exchange is made possible using just two channels uses, 
unlike the existing work done for the case, to the best of 
our knowledge. The Relay node clusters the 4^ possible 
transmitted tuples {xA^ xb,xc\xi)) into various clusters such 
that the exclusive law is satisfied and singular fade subspaces 



are removed. This necessary requirement of satisfying the 
exclusive law is shown to be the same as the clustering being 
represented by a 4-fold Latin Hyper-Cube of side M, if M-PSK 
is used at the end nodes. Using the proposed scheme that uses 
two channel uses, the size of the resulting constellation used 
by the relay node R in the BC phase is reduced from 4^ to lie 
between 64 to 90 for M=4. The size of the clustering utilizing 
modified clustering may not be the best that can be achieved, 
and it might be possible to fill the array with lesser than 90 
symbols. 
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APPENDIX 



Null Spaces of the Singular Fades Spaces for the case xa — x'^, xb — x'q, xc — x'q and xd — x'j-, GVi. 
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